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ABSTRACT 


азис Ог а spectral formulation of the primi- 
шише сашаттопз using a diagnostic divergence is tested for a 
global model. The initial conditions are generated from a 
developing baroclinically unstable wave. A semi-implicit time 
scheme is developed and tested along with the usual explicit 


method during the course of the experiments. Results show a 


лат ve ly small effect of a divergent initialization on the 
ensuing integrations. The semi-implicit method shows a ten- 
ШЕШЕ TO стостп OU High frequency oscillations in local ten- 
dencies. 
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Б ODUCTION 


The emphasis of this thesis is two-fold. The first point 
Пи СЗ осе ово ime characteristics of a spectral 
primitive equation model. The second is to attempt to Таке 
advantage of the spectral expansions in an initialization 
Technique. 

A spectral model uses a Galerkin formulation with the 
еса! harmonics as its choice of basis functions. The 
Basis functions are eigenfunctions of the Laplacian operator 
in spherical geometry and, thus, are used to advantage for 
global integrations. An immediate consequence of the Galerkin 
formulation is that spatial derivatives of a single wave are 
computed without the usual truncation error present in finite 
difference methods. This means that the spectral model 
contains no linear phase speed errors. 

The non-linear terms are computed via a transform grid 
in physical space (Orszag, 1971). This procedure is respon- 


sible for making the efficiency of the spectral model com- 


parable with finite difference methods. Mao Ee 
non-linear terms may be computed without aliasing. This pre- 
Bm ШЕЕНІ Іпетасі |ІТту апа conserves Invariant integral 


properties. Physical processes are also computed on the 
means 1 ORM arid. 
The second objective of this thesis is to test a divergent 


шша галоп Гог the spectral model. Most static and 
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— K 


ре Ii aiVzatren Techniques of global primitive 
equation models use a constraint to balance the rotational 
part of the wind with the mass field while setting the 
divergent part to zero. Phillips (1960) has shown that a 
completely non-divergent initialization cannot eliminate the 
gravity waves which would be generated and his work suggests 
that the quasi-geostrophic divergence could tend to suppress 
these gravity wave modes for mid-latitude synoptic scales. 
С Иса пе пле quasi~geostrophic divergence had interna! 
modes close to those which would be developed by the model, 
such things as the large scale precipitation would be more 
NAS Tic during The first few hours of the model forecast. 

u amor This Thesis TO develop a diagnostic 
velocity divergence to insert into the fields of an otherwise 
completely balanced system and to test the subsequent effect 
on The spurious gravity wave noise and the developing diver- 
gence. 

ee ral formulation of The primitive equations can 
SET CG echni ges which initialize according to scale. 
For example, the quasi-geostrophic divergence might be com- 
puted only for the synoptic scales of the model. Furthermore, 
the model conveniently supplies diagnostic information about 
the various scales. 

The final goal of this thesis is to develop a semi-impli- 
Cit time differencing scheme and test its effects during the 


marsi fenm OUr S after initialization. 
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The model 


NOBEL DESCRIPTION 


described here has been developed by Rosmond 


(1977) and the details are presented here for reference 


purposes. Similar formulations have been done by Hoskins 


MaS immons (19753) and by Bourke (1974). The equations for 


ЕШ іпуіссіс, 


written as 


се; 12 
HG +I 


е 
v 
Q2| Q2 
+D 
и 


is 
+ 
ч" 
Q| ap 
= | О 
tl 


where 


< 
I 


adiabatic and hydrostatic atmosphere may be 


> ~ · ду 
-Ve (C+f)V — Kk ° Vx CRTVq a с 55) 
> >? 
ke Vx (t+t)V zov Уб FT se - v (à N ) 
да 2 
> 90 
-V VO - 0 за 


= да 
-0 -у Уа 7 $0 


= Т 


vorticity ( = k * V x V) 
velocity divergence (D= V >œ V) 
Temperature 

potential temperature 

surface pressure 

horizontal velocity vector 
geopotential height 

gas constant 


specific heat 
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f - Coriolis parameter 


O - vertical coordinate (c = р/т) 
б - vertical velocity (0 = do/dt) 


q = пп 


The continuity equation (1.4) may be refined by integrat- 
ing in the vertical and imposing boundary conditions on а; 
ШІ) со) = 0. ШЕР еа Ире оона топ G = V*Vq, we may 


демо те Eq. (1.4) as 


9q _ - = 
б) ugs m (D + G) 


where the over bar denotes a vertical average. 

The vertical velocity, on may be obtained diagnostically 
ВЕ шта Eq. (1.6) into Eq. (1.4) and integrating in 
шише vertical to obtain 


G 
|.7) с- 6 + ба - / (G + D) da . 


O 


By defining a horizontally mean temperature, 
х 
Il CO) | (а. А, 9, t) 
ЕШ Мета The following operator 
oO ca 


a(X,Y) = a т 3u 
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We can now expand the vector quantities in our basic 


equations in spherical coordinates (A,n) to arrive at 





oC А 
ЕС) ЭТ = (А,В) 
* 
1.9) = = a(B,-A) Es Га) 
99 _ . 99 
0) дн 17 a(U8,V8) + 8D - С 55 
Jo = 
Ш) m = -(D + G) 
I2) с ELE -RT 
да 
мһеге 
> * ОМ RT' 
A = (C+f)U + б 55 + — cos Ө 
r 
к . 9) ВТ! 9дпт 
В c ат 3À 
и 
G = U дит да пт 
а дА дъ 
E u2 + ү2 
Е ne 
2(l-u ) 
U - u cos n/r 
MESE OOS 
А = jongiTude 
n = latitude 
u = sinn 
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аси CIZ) are The basic equations used in 
the model. The equations are represented spectrally in the 
лопта | and и ике ditterenced in the vertical. The de- 
pendent variables are written in terms of a triangularly 


truncated series of spherical harmonics: 


where the summation is a double sum over m,n for |m| < М апа 
Im| < £ « M. The separation is such that the coefficients, 


X ае типсттоп5 от Time and vertical coordinate and the 


m 
2? 


spherical harmonics, Y s oum |! рипгетјоп5 от space. 


m 
02 


ШЕ Погпа! ігдатіоп апа огтподопа 11 ту pmeperties of the үт 


И 


Бои the coefficients to be obtained as follows: 


2T + | 
Xp - «X, (1% Б Іт / / X v арал 
О = 
The model used for this study has five vertical layers 
Staggered as in Fig. |. The even levels carry о апа the odd 
levels carry all the dependent variables plus the diagnostic 


Vamtanles, U,V. 
The finite difference scheme for the vertical advection 


terms is 
OX | М ° А 
Я 90 . 2 (A9) LOG LI ДОО Кы a E EN k 


This formulation conserves kinetic energy. 
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( 


Following Arakawa and Lamb (1976), the hydrostatic equa- 


mun (1.12) can be finite differenced as 


(8 
" es 
БЕ) 70-0 а 010 
where 
A л R/C РА AM R/G 
а | -(0,/0, , |) р and 8, - (0, + I P 


E cumdarv condition is formed by integrating the hydrostatic 


Sota ¡on from c » O fo | which gives 


or 


where is The Terrain geopotential. This boundary con- 


Va 


gon is then finite differenced as 


(До) 


|. 14) с Ф (да) - ПЕЕ к Т, 


Ом (1.13) and (1.14) are then combined to give the 


inso equation 


AE A 


which is The finiTe difference form of the hydrostatic equa- 


Tlon where the variables are now column vectors. 
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The pressure tendency equation (1l.1ll) can be written as 


eO N 
LO, T D, (Ac 


) 
ot К K 


ШЕ Іп matrix form as 


ше ea = -[М] (6+0) 


“ООШ "еашаттоп (1.7) can be written in matrix form as 
ГІТ) o = [2] (G+D) 


The thermodynamic equation (I.!10) is differenced consis- 
tent with the method employed in the hydrostatic equation. 
The potential temperature is written in terms of temperature 


and the last term becomes 


G с R/C 
S | я k _ R/C . E ро 
a comete ke e PTS TC. - 
k O 
K+ | k-] 
RT, 2 2-2 
+ سج‎ )6G = G = (۰ 
р 
For The purposes of the semi-implicit formulation де- 


scribed later in the paper, we separate the temperature into 


x 
its horizontal average, T , and perturbation, T'. The 


Ds ЕТ! Із defined such that 








Е б, R/C E б, R/G Í 
T |C Y, Т а: TES 
kt | | = -——— 
Z (Ao) n — 2006) . o 7 L6 
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The vertically averaged divergence can be written D = 


MO, so that 


pe m па О т АТ мо. 
oe р 


 егітітіпд for o from equation (1.17) we obtain 


-[YICZJe -([YICZ] + 2 [NDD 
Ср 


ORE 
-LY Jo Ue LN.]D 


-5 16 + Г9 0. 


Finally the entire term can be written as 








Or R/C, 
e - ! e ! 
1.18) = -[sJo+[QJo- 7087 - [6, С С io 
G 
Kt | 
С, R/C, A 
б ! = = = 
ИЕ Ир со D 
о i P 


The nonlinear terms are computed using the Transform 
ШЕШЕ 5йсдевстесй by Orszag (1971). The longitudinal direc- 
Bes dene using Fast Fourier Transform and the latitudinal 
direction is done with Gaussian Quadrature. The number of 


des М, апа longitudes, M, satisfy: 
2 па М> 5111. 


The number of points are picked to ensure non-aliased results 
for quadratic terms. Certain terms containing Об have triple 
products and are not computed alias free. Experiments by 

Шо О О аы а previous study (Hoskins and Simmons, 1975) 
indicate that this source of error is negligible for the 


length of the integrations performed in this study. 
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Пи oer IMPLICIT SCHEME 


This section gives details of the semi-implicit time 
differencing scheme used in some of the experiments of This 
study. Using the matrices defined in the model description 


urn on,seauations (1.9)-(1.12) can be written 


ж 
ШІ = a + RT q) = F 


br D 


3T I 
2.2) Fe - 1910 = ғ; 
29 94 = Nie = -с 
2.4) Ф" = ГС 
! = 
where ф ф Я 


Ғы апа Ғ. represent the remaining terms in each equation 


D I 
which are not explicitly separated out. Following Robert 
а 19727), all the terms on the left hand side of equa- 
tions (2.1)-(2.4) are time averaged. We define the following 


averaging operator and represent the local tendencies by a 


Leapfrog finite difference scheme: 


ot ƏК ГАТ пера 41 
2 
БИ ET V. ( ‚trat - ( 1751 
dT t 2(At) 
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è 


Бисел С221)/-72.4) Can Then ре written as 


— хо. 
пе 6D = Е - е + ВТ a) 
if D 
тг о 
| + Т 
a = 
РК © [N] D 


ПИ  пеј T 


Eliminating m ПОЗ ОООО Ко У св ти топ of Eq. 


(2.8), we arrive at a time averaged set of equations: 


E - X 
2.0 HALF, = о =) 
Ten Han (E, + [QJ 51) 
ПИ = oi ONG + NID 


СЕЛИ па equations (2.11) and (2.10) into (2.9) and 
solving for D' we obtain a Helmholtz equation for the time 


averaged divergence 


t _ nt-At 


2.12) (СИНА led =D + E CN 


t-At 


ж X 
#САТ) [С ЈЕ ART 4 -(At) RT С) 


Eee cent, matrix ana [B] = CCICo] - RT CN]. 


oO e Troll ormulalion, Ihe equation iS Trivial to 
solve since the spherical harmonics are eigenfunctions of 


the Laplacian v^. Knowing gu we then compute the prognostic 


2 








variables po "upon lascueemguafyions (2.6) - (2.7). РАТ 
is computed from the definition of the Time averaging 
се гатог: 
pt tat £ 257 _ 11741 
ЕП Тһе vorticity, Ani, never entered the semi-implicit 


formalism and is computed from a standard Leapfrog scheme. 
Calculations with the semi-implicit model required 5$ more 
eempuier time per Time step. Stable integrations were 


possible with a time step of 60 mins. 


p 





RA E EOS TROBAR DIVERGENCE 


КО ОУ Ге analysis of Simmons and Hoskins (1976), a 
mumunostic divergence Ís obtained for the case of quasi- 
ле Fophic flow with the full variation of the Coriolis 


parameter over the sphere. The relevant equations become 


- V, * V(c*f) - LD 


br) Y 


аха» 
+“ 
|| 


x 
3.2) Xz - V (0! + RT q) = O 


x x 
ol mE р от 
22) ao DT D e s 55 ) 
мһеге Ya Ета расову 
5 - sigma dot at odd levels 
Lo s V (HOT ^D) 


Making use of the hydrostatic relationship, equation 


(1.15), we may write equation (3.3) as 


ao! _ _ Е Zu; 
5 = “1y ° тт) +£ y Š 
И x " 
where у = [С] E I - = д, 
р 


Taking 72 Ae lion кожног! equation (3.1) and $— 


weed. oe sand SUbSTITUTINg we arrive at the following 


A 


equation for 6: 


5 -2 2 рый Од 
м L (Vg УСЕЊЕ) )+ D. HEC JC, VT) RT VS 
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A second equation for G is obtained from the continuity 
equation in the form of eq. (1.17), which when explicitly 


тей out, is 


Е 5) с - а(0 % С) = + (р + б) да 


where с is specified on even levels. 

The diagnostic x from equation (3.4) can be used to com- 
pute the divergence using eq. (3.5) provided we specify a 
HEsugsrv condition. For the purposes of initialization, this 
is most conveniently chosen to be 


A 5 = 
Е = UD cp) 0 


ens will tend To eliminate external gravity waves. 


ШЕ (Па this Constraint to equations (3.4) and (3.5) we ob- 


тап 
£4 -2 
3,6) үбс--7 (Vy Gee 20) + [Cave 
A С 
3.7) С = -fo + G) do . 
о 
Тһе method of solution will be to compute O using equa- 


tion (5.6) from an initial guess of the divergence. The even 
level o's are then found by interpolation and a new guess of 
the divergence is obtained by inverting equation (3.7). The 
updated divergence can be substituted back into equation (3.6) 
and the procedure continues iteratively until sufficient 


accuracy as obtained. 
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Бер 220) and (2.7) will now be put into spectral 
form. The equations contain three similar advective terms. 
Ex be sufficient To show the spectra! form of one of 


Them: 


Ve. A 


5.8) Y 2 ЭШ ЭА DA Ou 


Y 


MS the Terms are individually transformed to physical 


space by the summations: 


M 
oW [a А 
пи Tu 
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ЗА m л 2 
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ды) _ m & ,imA 
==>) Wt tm. 
£,m 


ШЕ ти! Тр 1еаттопз in eq. (5.8) are Then computed at 
отт end The result is transformed back їс spectral 


space by computing the inner product defined as: 
+] 2T 
Ше 1 М 2 M 
529) Jy = = | | JO ANDY dunes rrt), U : 
-| 0 


The integral is computed using Gaussian quadrature in 
The Uu direction and Fast Fourier Transforms in the à direc- 
Is ТЕГІ Of equation (5.6) can then be 


evaluated: 


б .(#7(02)) 
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and expanded in spectral form using the previous computation 


m 
of Jo: 
И +1 D ee 
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where the à integration has already been done. 


Making use of the recurrence relationship 
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we can write the first integrand as 
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and using the orthogonality of the Legendre polynomials we 


Samen for the first integral 
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The second integral can be done similarly using the 


relationship 
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To obtain 
| М 


2 J ЭР 
ip f aa Ди 


> 
, 201) du 
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M M M 
мм Ян M EL 
ШЕ ТЕЗ ЛЕР 
=! 


ona ess olaaa SI) we obtain 


Ф T — — + таарч. 


The second order operator £D can now be formulated. 


Using The nofation, < ZD, E т 2i , we can write by analogy 
ШЕП equation (3.12) 
5а еМ pM L М M 
BN EU E 
Шише”, 
END ан 
Dr 
em 
and 
Z M _ M = a M Base M 
Ec О Se Ft er ey, ud 
я : : M M ; 
атола in for Би апа Play we obtain the result 
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КОШ e 310 (2.7) may now be fully expanded 


as 
<M | М a M 2. „М 
5.14) 6 + p * Eae VE) Y > + <£%D, Y >) 
M У М М 
ЕЙ ар = -/ (0 + 60) do 
О 
where 
M _ | М 
DE Be} an ШАП. 
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eevee Or) ine EXPERIMENTS 


The initial conditions used in the experiments of this 
thesis were developed from a baroclinically unstable mean 
flow and a small perturbation in zonal wave no. 6 (m=6). 

The mass field of the mean flow was constrained to be in 
geostrophic thermal balance with the mean wind. The vertical 
profile of the mean wind was linear varying from 4.5 m/sec 
ат с = | to 49.5 m/sec at o = 0. The latitudinal structure 
Varied as sin’(2n) placing a jet max at n = 45°. The verti- 
Бай етоегатиге рготі|е approximated the U. S. Standard 
Atmosphere. 

The perturbation in wave no. 6 corresponded То а maximum 
v-component of „| m/sec. The model was integrated for 192 
hrs allowing the baroclinic wave to grow to finite amplitude. 


Er shows a time series of the vorticity during this 


Ai The growth is strikingly linear until 192 hrs where 
there is evidence of the mean flow (m = 0) changing due to 
feedback from the eddies. The linear growth corresponds to 


md time ot raughly 26 hrs form = 6. The fields at 
192 hrs provided the balanced fields for the initialization 
experiments. 

Three types of experiments and one control case were 


integrated for 72 hrs of model time. 
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и 


Пт нет сео =е consisted Of initializing with The 
balanced fields and integrating. The experiments consisted 
Еш The following: 


Experiment A: Initialize wi th The balanced fields but 
no divergence 


Experiment B: Пипина те ми" the balanced fields plus 
the quasi-geostrophic divergence 


Eel Ment C: ШЕҮп езі to Experiment A but with a semi- 
implicit scheme 


All experiments and the control were first integrated 
МІІП а spectral truncation of wave no. 2l. These will be 
referred to as low resolution experiments. High resolution 
experiments were integrated for the control and Experiments 
А апа В. These cases were truncated at wave по. 42. Time 
Sieps of |2 mins for the low resolution and 6 mins for the 
high were used. Experiment С was integrated with time steps 
of I2 and 60 mins. All experiments had a weak Robert time 


Mero 0. | and no dissipation Term was used. 


а 





A RESULTS 


CONTROL CASE 

s c 3-6 Show the initial fields at level 5 for the 
low resolution experiments. Figures 7-10 are the correspond- 
ing high resolution fields. They are characterized by a 
ВЕ орпа frontal zone. The advantage of the high resolu- 
tion is apparent in the temperature gradients across the 
ШИ ГОП апа warm front amd in the surface pressure. The 
low in the high resolution surface pressure field is 5 mb 
deeper. The negative and positive areas of fhe divergence 
and vorticity fields are approximately symmetrical in the 
resolu ont 1elas and quite asymmetric in the high reso- 


Цикоте јаз. The region of maximum positive vorticity 


corresponds very well with Тһе maximum temperature gradient. 
The mean fields consist of an indirect cell induced by 
B -adgrowing baroclinic eddies. The indirect cell develops 


In order to decrease the vertical wind shear which must remain 
in approximate thermal balance with the decreasing pole to 
ERN temperature gradient. The effect is more pronounced 
in the low resolution fields. The nigh resolution mean 
Misco nta in a чесестпаггест сігсиіатіоп апа а strong 
тест cell. 

iene monte sittarion during The period of the control 


ы СШ СЫ Ое О а rapid cascade of energy to smaller scales. 
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СЛО ОПТ 5 сафе Of The frontal zone decreases. Figures 
li-15 show this effect in the high resolution temperature 
needs for Times out TO 56 hr. Beyond this point, the 

energy builds up at The spectral limit of the model since 
mene is no dissipation and the 72 hr fields are non-meteoro- 
logical (Fig. 14 and 15). Because the energy blockage occurs 
al smaller scales, the high resolution case (Fig. 14) is, 


paradoxically, noisier and in greater error. 


Bee GRAVITY WAVES 

Мат пе а оп of later sections, we discuss the 
Characteristics of inertial gravity waves for the model used 
ШІ ТІПІс study. The E e rods is The most sensitive indica- 
Tor of gravity wave motion. As discussed in Hoskins and 
ое), equation (2.12) 15 тһе semi-implicit analogue 
of the gravity wave equation for a multi-layer model. The 
eigenvalues of the matrix [8] (eq. 2.12) give the permissible 
eigenvelocities of the gravity wave modes. Тһе corresponding 
eigenvectors give the vertical structure. Figure 16 depicts 
The eigenmodes for the experiments in this study. Table | 
gives the periods of the gravity wave modes for zonal wave 
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шосета но е сепсе тау be Thought of as consisting of 


a sum of the ¡individual eigenmodes: 


4.1) и-пи а В 
k 


where the D, are the eigenvectors. Тһе weighting coefficients, 


K 
a,» may be obtained from the simultaneous equations formed 
Dune dot products of the individual eigenvectors with equa- 
wen Al). A Time series of the coefficients will provide 
ШЕГІП! dragnostic information іп |атег sections. Figures |17 
EID are da time series of the external and first internal 
modes for the low resolution control case at a point on the 
at itude circle. It is apparent that the coefficients 
contain both the meteorological modes and any excited gravity 
memes. Figure |/ contains no high frequency oscillation and 
is purely meteorological. Figure |8 contains a meteorologi- 
оде мт а period of roughly 5 days plus a smaller 
amplitude wave whose period of approximately 1/2 дау согге- 
sponds well with the calculated value for the first internal 
mode from Table |. The remaining internal gravity modes 


have periods approaching meteorological values and are diffi- 


ӘШІР 10 distinguish. 


OE RERI MENTS A AND B 
low Resolution 
Neglecting the divergence in the initialized fields 
(Experiment A) generated computational gravity waves.  Compar- 


ing the pressure tendencies of the control case (Fig. |9с) 
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ose об Ехрепштепт А іп ГІ. |9а, ІТ іс apparent that 


ШЕШІ (гесцепсу Oscillations Of roughly 5 hr period are 
excited. The ¡initial magnitudes of the oscillations аге 
approximately 1-2 mb/hr. The oscillations dampen out with 


time as the geostrophic adjustment takes place although the 

dispersal of the gravity waves during the adjustment process 
is hampered by the longitudinal periodicity and the symmetry 
of the fields about the equator. 

Figure 196 shows the pressure tendencies with the 
quasi-geostrophic divergence inserted into the initializa- 
tion fields (Experiment B). The highest frequency oscilla- 
tions have been reduced over. Experiment A by a factor of 
approximately 2. 

Comparison of The divergence modes in Figs 20, 21, 22 
allows us to be more quantitative. The figures have the 
Sommmo! case subtracted out and contain the gravity wave 
noise generated by the approximations of each experiment. 
The noise generated in the internal modes by the non-divergent 
initialization is an order of magnitude larger than the ex- 
ternal mode and is generally 25% of the total divergence. 
The figures show that the external mode of Experiment B is 
about 68% of сураи О. А. The first internal mode shows 
more noise in Experiment B initially and then decreasing to 
50% and the second internal mode shows a decrease of about 
50% over Experiment A. In general the effect of the quasi- 


geostrophic divergence is to reduce the gravity wave noise. 
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25 Птап Resolution 

Figure 23 shows the quasi-geostrophic divergence at 
level 5 which was used for Experiment B. A qualitative com- 
parison with the actual model divergence (Fig. 10) shows the 
phases to be in agreement but that The quasi-geostrophic 
approximation is an overestimate for the perturbation fields. 
The zonal mean fields do not compare well. The quasi- 

Пеште moph!ic divergence shows a strong indirect circulation 
which is displaced poleward relative to the model divergence. 
There is no evidence of a direct cell. 

Figures 24-26 depict the noise generated in the 
divergence modes for Experiments A and B. Generally they 
show that the gravity wave noise of Experiment B is 66%-72% 
Anar in ner: A. The effect of the quasi-geostrophic 
divergence is to reduce the noise. However, this reduction 
15 less than іп the low resolution experiments. 

In view of the fact that the high resolution experi- 
ments were less successful than the low resolution and that 
quasi-geostrophic formulations are valid only for small 
values of the Rossby number, it is interesting to compare the 
initial model divergence of the control case with the quasi- 
geostrophic divergence in a spectral manner. Referring to 
Fig. 27, it is seen that there is good agreement between the 
model and quasi-geostrophic divergences for zonal waves 
corresponding to m = 6, 12, 18. The overestimation by the 
quasi-geostrophic divergence noted earlier is primarily in 
wave no. 6. The mean fields (m=0) and waves corresponding 


ШЕП 24 О 56 42 do nof compare well. 
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One motive for initializing with the quasi-geostro- 
phic divergence is that, presumably, the "true" divergence 
ШЕШІП ре appreximated ar fhe start and that the model would 
mms develop This divergence field earlier than a non- 

En crgent initialization. Figures 28 and 29 address this 
zu They show that at 6 hours into the integrations, zonal 
ле штрег5= б, 12, [8 do, in fact, approach the control 

Ки гаепсе slightly faster. There is almost no difference 

for the higher wave numbers and the mean field shows poor 
adjustment in either case. This might be expected for the 
mean field since the geostrophic adjustment process ¡s very 
ШЕШШ степт ат this scale їп а bounded region. Оп the other 
hand, the higher wave numbers adjust very rapidly and after 

6 hrs, they are in good agreement with the control case. 

In light of the results showing good agreement in 
и WiWal quasi=geestrophic divergence with the control 
divergence for zonal wave number 6, 12, 18 and bad agreement 
Toe Scales, if IS natural to modify Experiment B such 
ШІСІ ӘПІу scales corresponding to m = 6, 12, |8 are used 
in The quasi-geostrophic divergence. Figure 30 shows Тпе 
time series of the external divergence mode for this experi- 
ment. Figure 3l shows the spectral adjustment of the diver- 
Qs а! © hours. The corresponding results for the unmodi- 
fied Experiment B are Figs. 24 and 29, The results are 
almost identical to the unmodified quasi-geostrophic diver- 


gence results. 
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Figure 30. Time series of the gravity wave noise 


іп Тһе modified Experiment В. Vertical 
scale is sec”"!, Time is un mr s 
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Be EXPERIMENT C 

ООО еони Ss Rias and Simmons (1975), a linear sta- 
Fury analysis Of the semi-impliciT scheme shows that The 
time step is limited only by the Rossby waves. The effect 
of the scheme is to slow down The faster moving gravity wave 
modes so that the integrations do not amplify in Time. 
Slower moving atmospheric motions are less affected. Table || 
presents sample calculations of the distortion of the periods 
of the gravity waves present in these experiments for two 
Time steps. These may be compared with Table I. The general 
behavior is that the periods increase with increasing time 


steps. 








TABELCE SHTI 
ое атта ашы ша 
ws mins Ж. p см 0 





Two points will be investigated in this experiment. The 
first is that the semi-implicit scheme should tend to smooth 
zh trequency oscillations, particularly with a large 
time step. The second is that the semi-implicit method 


might possibly hamper the geostrophic adjustment process 
since it slows down the gravity waves which are the mechanism 


for the adjustment. 
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Figure 32 shows the pressure tendencies using the semi- 
implicit method with time steps of 12 апда 60 mins.  Compar- 
Mais wich the corresponding results for the explicit 
ааган ол5= (Fig: 19a), if is evident that the effect of 
memsemi-implicit method iS to smooth out the 5 hr oscilla- 
ШИ Figure 33 shows the relative effect of a small vs. 
large time step in the semi-implicit method on the external 
gravity waves. The larger period of the 60 min time step is 
evident. The larger time step initially amplifies the com- 
рими топа noise, Бот а тег 24 hrs it is considerably damped. 
Тһе smoothing in the pressure tendencies is postulated То 
come about from the decreased velocity of the external 
gravity mode assuming an advective Time scale. 

The second point of this experiment is shown in Figs. 34 
and 35. Comparison of the divergence coefficients for time 
steps of 12 and 60 mins shows no difference at 3 hrs. The 
adjustment of the divergence to the control case ¡s apparent- 


ly the same for either time step. 
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Vil. SUMMARY AND CONCLUSIONS 


Китен тези Што ог" тако study are based upon rhe divergence 
developed by a baroclinic wave on a uniform earth. Dans ali = 
zation without the divergent wind component generated spurious 
inertial gravity waves. The internal state of the initial 
fields had nearly zero mean divergence. Thus, The externa! 
mode of the generated gravity waves was an order of magni- 
ae less Than The Total divergence and was not significant. 
However, the internal gravity wave modes generated were signi- 
ficant. An important source for the external mode, namely 
mepography, was neglected. | 

Initializing with the quasi-geostrophic divergence 
generally decreased the spurious gravity wave noise and aided 
the ensuing integrations to develop a divergence closer to 
the control case, but both effects were small. Furthermore, 
the advantage of the quasi-geostrophic divergence decreased 
when smaller scales were included in the Integrations. In- 
serting the quasi-geostrophic divergence into only synoptic 
scales gave no advantage. However, the effect of this experi- 
ment may have been reduced by the fact that there were no 
planetary scale waves. 

Using a semi-implicit scheme with the same time step as 
an explicit integration had no effect on the divergence modes 


but did smooth out oscillations in the pressure tendencies. 


T3 





A larger time step initially increased the amplitude and 
вто ог she fast external gravity waves and had less effect 
on slower moving mofions while still stabilizing pressure 
tendencies. 

In general, the effect of a divergent initialization for 
a global primitive equation model is small. This same con- 
clusion has been reached by two similar studies = Houghton 
сие (1971) and Dey ef al (1975). This study, however, has 
neglected any physical processes such as precipitation as 
СЕЛЕ А5 Orography and the effect of an initial divergence in 


These cases may well be more important. 


74 





RE ERENCES 


кома. A., and V. R. Lamb, 1976: Computational Design of 
the Basic Dynamical Processes of the UCLA General 
ee ularion Model, 160 pp. 


Eme W., 19743 A Multi-Level Spectral Model I. Formula- 
ШІСІ апа Hemispheric Integrations. Mon. Wea. Rev., 102, 
687-701. 


EIN а. A. Brown, and R. D. McPherson, 19/5: Some 
Initialization Experiments with the NMC Global Model. 
The GARP Programme on Numerical Experimentation, Rpt. 
ie. 8. 


по, В. J. and A. J. Simmons, 1975: A Multi-Layer 
ес га мее аа The Semi-Implicit Method. Quart. J. 
а METE. SOC, 101, 657-655. 


Вог, 05 0. Бапппетпег, and М, М, Washington, 1971: On 
Mota izarion of The Primitive Equations: 
Part 11. The Divergent Component of the Horizontal 


пирон | Appl. Метеог., 10, 626-654. 


MS. Ac) |971: Numerical Simulation of Incompressible 


Flows within Simple Boundaries - Galerkin Spectral 
EEDuesensatioms. Studies in Appl. Math., Vol. L, 295- 
27, 


Sd.) J. Henderson and C. Turnbull, 1972: An 
IA (Іле Integration Scheme for Baroclinic Modes of 
the Atmosphere. Mon. Wea. Rev., 100, 329-335. 


IN anda В. 4. Hoskins, 1976: Baroclinic Instability 


on the Spheres: Normal Modes of the Primitive and Quasi- 
БЕРЕБШЕСЕШІС ЕШшатіопе. J. Atmos. Sci., 33, 1454-1477. 


75 





L HE sns DISTRIBUTION 


Defense Documentation Center 
Cameron Station 
Alexandria, Virginia 22514 


Library, Code 0142 
Naval Postgraduate School 
Monterey, California 93940 


Мо | Halmımer, Code 65Ha 


Chairman, Department of Meteorology 


Naval Postgraduate School 
Monterey, California 93940 


ШІ 5. T. Williams, Code 63Wu 
Department of Meteorology 
Naval Postgraduate School 
Monterey, California 93940 


тесто 
Naval Oceanography and Meteorology 


" 


National Space Technology Laboratories 


ES ours, MIssisslppi- 59520 


Siricer in Charge 


Navy Environmental Prediction Research 


Bacal Ty 
Monterey, California 93940 


Commanding Officer 
Fleet Numerical Weather Central 
Monterey, California 93940 


Sas |І, Rosmond 


Naval Environmental Prediction Research 


Facility 
Monterey, California 93940 


mu molut M. Lubeck 


COMNAVMARIANAS Box I2 
нео ап Francisco 96630 


76 


Мо. 


Соріев 


N 














MB 1 o mE 








тын = марла 
ia en DE | 
ME 
x E u I | B 
I = 
Iz С u I mE NE Bo dm Tos > — I 





Т 


и 

















Thesis 4 74) 
1885 Lubeck . 1188! 
ca Divergent initiali- 

zation experiments us- 

ing a global primitive 

equation spectral model. 








